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Abstract (In j. 


An identification procedure is proposed to identify the damage characteristics (location and 
size of the damage) from dynamic measurements. This procedure was based on minimization of 
the ’mean-square’ measure of difference between measurement data (natural frequencies and mode 
shapes) and the corresponding predictions obtained from the computational model. The procedure 
is tested for simulated damage in the form of stiffness changes in a simple fixed-free spring-mass 
system and symmetric cracks in a simply-supported Bernoulli-Eulcr beam. It is shown that when all 
the mode information were used in the identification procedure it is possible to uniquely determine 
the damage properties. Without knowing the complete set of modal information, a restricted region 
in the initial data space has been found for realistic and convergent solution from the identification 
process. 


Introduction 


There is a considerable body of research on identification problems, that is, the problem of 
identifying the engineering properties or reconstructing the structural configuration of a vibrating 
system from certain natural frequency spectra and/or corresponding mode shape. Such problems 
were considered by Barcilon [1, 2], McLaughlin [3, 4], Gladwell [o-<], and Gladwell ct cil. [8J. Most 
of these studies involve the determination of material properties from natural frequencies, and 
they emphasize the existence, uniqueness, and methods for determination of properties (tcimcd 
’reconstruction’). 

An detection procedure was developed by Shen and Taylor [9] to determine the crack character- 
istics (location xc and size cr of the crack) of Bernoulli-Euler beams from their dynamic response. 
The idea of this procedure was related to methods of structural optimization. Specifically , the 
structural damage was identified in a way to minimize one or another measure of the difference 
between a set of data (measurements) T d , and the corresponding values for dynamic response M d 
obtained by analysis of a model for the damaged beam. This may bo expressed symbolically as the 
following optimization problem: 

min norm(Td - M d )- ' > 

xc,cr 

Naturally, the minimization represented here is constrained by the equations which model the 
physical system. Moreover, as indicated in the discussion by Shen and Pierre [10, 11], one can 
note that the more modal information used for crack detection, the more accurate and reliable the 
result that can be achieved. For practical purposes, the objective of Eq. (1) was formulated based 
on a certain set of specific modes; specifically the first three modes are considered in the inverse 

procedure. . .. 

In this study the corresponding to the mean-square measure of the norm, as shown in Lq. ( i), 
is examined. The identification process is based on minimization of the ’mean-square’ measure of 
difference between measurement data (natural frequencies and mode shapes) and the corresponding 
[predictions obtained from the computational model. The identification procedure is tested forj 
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Emulated damage in the form of a symmetric cracks in a simply-supported Bernoulli-Eulcr beam 
and a fixed-free spring-mass model. The uniqueness and reliability of the identification process is 
confirmed by solving several damage identification examples with specified damage positions. 

Problem Statement 

In this section, variational formulations for the identification of damaged one-dimensional 
structures are presented. The mean square differences between measured and modeled values of 
frequency and mode shape are employed as the objective function in one of the formulations. In 
other words, the inverse process seeks to determine the damage parameters, location xc and size 
cr, in the mathematical model to minimize the mean square difference between the test data and 
analytical predictions. The problem formulations are presented in forms of a cracked Bernoulli- 
Euler beam and a multi degrees of freedom (DOF) spring-mass system. 

Cracked beam model 

In the treatment of this problem, it is assumed that the testing information (data) is provided 
from certain test points distributed over the structure. This data is comprised of frequency and 
mode shape information associated with the lower several response inodes. 

For a simply-supported uniform beam containing one pair of symmetric cracks (see Fig. 1), 
the problem of optimization in crack detection can be expressed, in terms of comparisons between 
modeled response and test data, as 

min [normal - w to {x tm ) - w a {x lm ))] (2) 

cr t xc 

subject to constraints that define the beam response iv a (ic., the equations for free vibration), and 
which prescribe appropriate normalization of w Q and test data w tCi . 

Here cr = —p represents crack ratio (a measure of crack depth), and xc identifies crack 
position (see Fig. 1). Also, the objective function measure of differences between measured and 
modeled values of deflection and frequency in Eq. (2) is stated for present purposes in the form: 

M T i 

norm{Jl Q - U 2 a , W ta - W Q ) = (^[(w^ - + £](u>la(Ztm) - Wa(*t m )) 2 ]) 5 (3) 

a=l m=l 

where w a , w a represent the natural frequency and mode shape of ath bending free vibration 
mode, M is the number of modes for which test information is available, and, once again, the 
corresponding test data are symbolized by and w tQ . Here x im (m = 1,2, ...,T) locates the m-th 
out of T measure stations, respectively. The measures w tQ and w Q that appear in the norm must 
be normalized on a common basis in order to facilitate comparison between the data and model 
values. 

The symbol $ is introduced to represent the square of the norm given in Eq. (3). The 
identification problem now can be stated: 


min $ (4) 

cr,xc 

subject to: 

f {EIQ(w'^x)) 2 - u>lpA wl(x)}dx = 0 (5) 

Jo 

T - 1 

^2(w a {xt m )v}p{x tm ))Ax tm - T ) al 3 = 0 

m= 2 


( 6 ) 
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P (cr + axe) - JZ < 0 (7) 

cr < cr < cf (8) 

xc < xc < xc (9) 

where a, /3=1,...,A/, a is a weighting factor on the cr and xc , /Z represents the upper bound on 
value cr + axe, and xc,xc, and cr,cr represent the upper and lower bounds of tire crack (damage) 
parameters xc and cr, respectively. (Note that both upper and lower bounds on the variables cr 
and xc are necessary in the present problem.) Since w Q comprise an orthonormal set, r) a p is defined 

aS T 1 

Waa = X) W l km Ax tm, li™ T ] a 0 = 0 for Q # P (10) 

1 1 T—*oo 

m = 2 

The effect of cracks on the structural properties of the beam is reflected by factor Q in Eq. (5), 
as described for symmetric surface cracks in Shen and Pierre [10]. In other words, the optimization 
parameters xc and cr cited in Eq. (4) enter the problem via Q. 

According to the K-K-T (Kurash-Kuhn-Tucker) necessary conditions for the optimization 
problem Eqs. (4-9), there exist Lagrange multipliers A Q , , A Q/ g, and T k which satisfy the following 
equations (the notation ’I*’ refers to solution points): 

A 0 > 0 


o 

A 

0 

< 


Ti [(cr + axe) — R] |»= 0 

(ID 

[r 2 (cr - cr)] |.= 0 

(12) 

[r 3 (xc - xc)] |.= 0 

(13) 

[r 4 (xc - xc)] |.= 0 

(14) 

[r 5 (cr - cf)] |.= 0 

(15) 

The solution must satisfy the following three equations as well: 


[2(u4 -«„) + \ a pAC 2 0 ] |.= 0 

(16) 

[(EIQw'^(x))" - ulpAw a (x)] |.= 0 ; x tm < x < x <(m +i) 

(17) 


X- 1 ot-l 

{— 2(tn* a (x) - w a (z)) + [2A acr tc a (x) + ^ A Q pwp{x) 
m = 2 P— 1 

M 

+ A a pW 0 (x)]Ax tm + 2 \ a [(EIQw'^(x))" — LJ 2 a pAw a {x)]} | x=X(m |,= 0 (18) 

0=a + 1 

Note that the above equation of motion (Eq. 17) is valid interval by interval over the span of the 
structure. 

Finally, the conditions for stationarity of $ w.r.t. the optimization variables cr and xc (ie., 
the optimality conditions) are: 

E X ° {EI J‘ |^K'(*)) 2 ^) + ri - r 2 + r 5 ] |.= o (19) 


| 

i 
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M r l 

[£ A a (Ei J g«(x)) 2 rfx) + r,fl + r 3 - r 4 ] |.= 0 ( 20 ) 

a— 1 

• The problem formulation for the numerical method-mean square criterion 
The purpose in this subsection is to re-state the inverse cracked beam problem with mean 
square criterion, Eqs. (4-9), in the following form that is more convenient for computational pur- 
poses. With the introduction of symbols f and T for convenience, the statement becomes: 


subject to : 


M T 

min D« ta £cr ) 2 + - 1t'a(2-(m)) 2 ] 

—1 a = 1 m=l 


T T 

[ a *Q E (^afatm)) 2 - Co 51 ( u ’o( 1 lm)) 2 ]Al| m - T a = 0 

m— 1 m=l 

T-l 

^ ^ (^q )) — T)af3 = 0 

m — 2 

0 < cr < 1.0 


0 < xc < 1.0 

where a, (3 = 1 V ..,A/, variable vector = {cr, xc, £ 0 , u; a (2- im )}, 


( 21 ) 

( 22 ) 

(23) 

(24) 

(25) 


El IT 4 


(2G) 


and 

T T 

To = [a 4 Q ^(WtaiXtm )) 2 - Za )) 2 ]A:E, m ( 27 ) 

m=l m=l 

Spring- mass model 

The spring-mass model to which the present identification procedure is applied is shown in 
Fig. 2. It consists of 3 masses connected by linear springs of stiffness defined by 


ki 


= fc( 1.0- 


<ft”» \3 

O > 


(28) 


where dm, is defined as a damage parameter at i-th spring. If dm t is interpreted to represent 
the same physical meaning as cr does in the cracked beam model, the system’s damage condition 
may be introduced by specifying a certain value to ’damage parameters’. For instance, according 
to Eq. (28), a damaged condition can be constructed in which stiffness drops 25% and 50% at 
the spring 2 and 3. This is accomplished by assigning the values dm 2 and dm3 to be 0.2743 
and 0.6189, respectively. In a sense, the spring-mass model can be viewed as a simple simulation 
analogy of the cracked beam, ie., both extent and location of damage can be represented in the 
model. The fundamental frequencies u>,* of axial vibrations are related to the mode shapes u, = 
(ui,ii2,t/3)f q = 1,2,3 through the equations: 


uf[K a ]ui - Ziiij [/)u, = 0 \i= 1,2,3 


(29) 


where Therefore, the present damage identifi 


ition problem can be stated as 


min 


3 


3 


]0(&« - c) 2 ) + 

.=1 j = 1 


(30) 
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to: 

iiJ[K a ]ui -tiufui = 0 ;»= 1,2,3 

uJuj-Sij^ 0 ; t,j = 1,2,3 (32) 

0 < drrii < 1.0 (33) 

Numerical Analysis 

The numerical optimization technique set forth in this study for vibrating cracked beam iden- 
tification problems is accomplished using the VMCON optimization package program (this imple- 
ments a sequential quadratic programming method). The VMCON program uses Powell’s algorithm 
which is an iterative scheme designed to converge to a point that satisfies the necessary conditions. 
Additional information regarding to VMCON is available in Ref. [12]. 

Cracked beam model 

The cracked beam model to which the identification procedure is applied is shown in Fig. 1. It 
is a simply supported beam of length / equal to 18.11 of it’s thickness 2c/, with uniform rectangular 
cross-section area A , and a pair of symmetric cracks of cr = 0.5 located at mid-span (xc = 0.5). 

Unless otherwise stated, the damage properties (cr and xc) of the simply supported cracked 
beams are identified by direct solution of the optimalization problems described in the previous 
section. The sensitivity to chosen values for the initial crack position xc are discussed later in this 
section. 

• Examples with position of the crack (damage) specified 

Consider the first example for crack identification, the simply supported cracked beam, for 
which the crack position xc is known. In other words, only the crack ratio cr is to be identified; there- 
fore, the variables in this problem are cr, £',s, and mode shapes w Q (x) (xj = {cr,£ a , tc a (x* m )}, x 2 = 
{cr,( Q , flert})- This simplified example problem with the crack position specified (xc = 0.5) is 
presented to demonstrate the concept of the crack identification procedure described in the last 
section. 

In this example, it is assumed that the dynamic measurements are collected at 9 test positions 
(T = 9) equally spaced over the span. The first and last test stations are located at the left and 
right supported end, respectively. Hence, the length of each test span Ax fm ,m = 1,...,T — 1 is 
determined to be . In structural dynamic testing, ordinarily only a relatively small subset 

of the theoretically available eigenvalues and eigenvectors can be measured accurately, ie., realistic 
information on higher modes is difficult to obtain from the measurements at a limited set of test 
stations. Only information from the first three modes is to be used as test data in the present 
identification process. Furthermore, according to the observations in Shen and Pierre [10], the even 
modes of a simply supported beam are not sensitive to a mid-span crack; therefore, in effect only 
first and third mode (a = 1,3) information is used to represent crack damage. 

Once again, the crack identification problem presented by Eqs. 21-25 is solved here with a 
specified value xc = 0.5. For given initial values of x, this optimization problem is solved to 
minimize the criterion F. The results of the cases with various initial conditions are shown in 
Table 1. In order to clearly compare the results, only the first three variables, £i,f 3 iCr, of variable 
vector Xj are listed in the Table 1. 

In Table 1, the top row denotes the assumed crack ratio and corresponding first and third 
eigenfrequencies. The symbol * denotes the expected optimal solution through the identification 
process. The first two column entries, indicate the fundamental and the third frequencies 

corresponding to the initial crack ratio cr which is given in the next column. The last three columns 
give the final values corresponding to previous entry values. These final values are obtained at the 
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stage where computation is terminated when the further optimal search obtains improvements for 
criterion F less than the specified tolerance (10£ - 5 Was adopted in the present study). Recall 
that for an uncracked beam cr is identically zero. Therefore, in this example, it is decided to start 
with the case of the initial value cr = 0.0 and for each case thereafter the cr value is increased by 

0 . 1 . 

From the results presented in the first case of Table 1, one secs that the parameters £i,£ 3 , an d 
cr were identified to be 0.84684, 70.1348, and 0.50033 from 1.0, 81.0, and 0.0, respectively. The 
mean square critera F was cut down from 118.13502 to 0.42440E-5. The maximum error is less 
than 0.5% of the test data for these parameters. The results are also quite impressive for mode 
shapes. In order to observe the global variance clearly, the initial, final, and testing mode shapes 
are plotted in Fig. 3. Three curves appear on each plot: the initial mode shape, the finial mode 
shape, and the mode shape from test response. The final mode shape on these plots agrees well 
with the test mode shape. This is expected and verified the accuracy observed from the results in 
Table 1. It can be clearly seen that accuracy of the mode shapes will worsen if higher mode results 
are to be predicted. Improvement can be obtained by an appropriate adjustment of the location of 
these test stations. However, a sensitivity analysis of the test stations with respect to the accuracy 
of the dynamic measurements is required. This is not considered further in the present study. 

In Table 1, rows 5 to 11 present the results for cases with initial cr = 0.1 to 0.8. The corre- 
sponding final point values listed in the columns 4-6 show that these cases exhibit, as expected, 
similar solution characteristics and accuracy. This provides a physical understanding of the geome- 
try of the solution set: for the inverse cracked beam problem with specified crack position, the mean 
square criterion of Eq. (21) is a convex function and it is bounded by the constraints of Eqs. (22- 
25). Hence, one may conclude that the convergence of the present optimization problem is obtained 
independent of the initial data chosen. In other words, as long as the initial data is selected within 
the problem’s feasible domain, an accurate and unique solution through the identification process 
is expected. 

Clearly the prediction of mode 3 shape shown in Fig. 3 fails to reproduce the expected sin 
curve. This is because the 3rd mode shape was plotted based on the deflections of the mode shape 
measured at only 9 test stations. While this reflects a limitation on how well mode shapes are 
portrayed, the quality of the final result for the identification problem is unaffected. 

• Simultaneous identification of crack position and depth 

The second numerical example deals with the crack identification of a simply supported cracked 
beam with unknown crack ratio and with crack position unknown. In this treatment, the variables 
in the optimization problem axe cr,xc, and mode shapes tv Q {x) (xj = {cr,xc,f a , ic a (x< m )}). 
Due to the limitations of the VMCON program, the examples that concerning with the testing 
mode shapes w t provided in the form of continuous functions are not shown in this subsection. 

The formulation of the crack identification problem (Eqs. 21-25) is tested again with both 
crack position and depth are assumed unknown. In the first few cases, the simulated dynamic test 
measurements are assumed to be collected at 9 equally spaced test stations ( T = 9). The first 
and last test positions are located at the left and right supported end, respectively. This example 
will be solved a second time using an increased number of test stations, to provide information on 
sensitivity of the procedure to the amount of test data. 

In Table 2, the top row denotes the assumed crack ratio, crack position, and corresponding 
first and third eigenfrequencies. The symbol * denotes the expected optimal solution through the 
identification process. The first column entry T denotes the number of test stations used to collect 
dynamic measurements. The second and third column entries, indicate the fundamental and 

the third frequencies corresponding to the initial crack ratio cr and crack position xc, which are 
given in the next two columns. The last four columns provide the final values corresponding to 
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[the previous entry values. These final values are obtained at the stage where the computation is 
terminated, when the optimal search obtains step-wise improvements of F less than a specified 
tolerance {\QE - 5 in the present study). 

Table 2 shows that cases with T = 9 have the final values of f close to <f*, but almost all of 
these cases have unacceptable final estimates of xc and cr. For instance, if the initial position 
is selected as xc — 0.4 and cr = 0.4, the values of xc and cr at the final iteration are 0.99789 
and 0.36289 which are approximately 98% and 28% different than the given test data. In other 
words, evidently the configuration with xc = 0.99789 and cr = 0.36289 is able to provide another 
minimum value of the criterion (besides the one associated with the expected result). This cracked 
beam configuration is shown in the solid curve of Fig. 4 . The mis-match between final and test 
mode shapes can be clearly seen. This observation confirmed the unacceptable error previously 
obtained in the comparison of xc and cr between the final and test data. Except for the case with 
initial cr = 0.4 and xc = 0.48 which provides less than 1% estimation error, the rest of the cases 
in Table 2 with 9 test stations are also found to have similarly large estimation error. Therefore a 
dependable solution in crack identification is almost impossible to achieve on the basis of the 9 test 
stations simulated measurement information on first and third mode response. This confirmed the 
observations in Shen and Pierre [10, 11], ie., for a cracked beam with an unknown crack position, 
a unique solution is not to be expected. 

However, by comparing the third mode shape in Figs. 3(b) and 4(b) to the mode shape in 
Fig. 11(c) of Ref. [10], it can be seen that an accurate third mode shape can not be approximated 
based on the displacements collected from 9 test stations only. This implies that the accuracy of 
the above computational identification might be improved if the third mode is approximated well. 
Therefore, the cases with more test stations should be examined since they would clearly provide 
better mode shape approximation. The largest number of test stations which can be accommodated 
in the identification procedure is 45, due to the limitations of the optimization program package. 
Once again, the test measurement points are equally spaced, and first and last stations are set 
located at the left and right supported end, respectively. The VMCON problem formulation is 
identical to the case of T=9; however, the variable vector x is expanded from 22 components to 94. 

Rows 12 to 17 of Table 2 summarizes the results through the minimization process. As in the 
previous cases, the final values of frequency £ are observed to be close to test values Acceptable 
final solution values for xc and cr are shown in the results of the cases in which initial xc and cr are 
selected within the range from xc — 0.4, cr = 0.4 to xc = 0.6 , cr = 0.6. On the other hand, within 
this range, good agreement is also shown in mode shapes. Figures. 5 and 6 display the initial, 
final, and test mode shapes for cases with the initial xc = 0.4, cr = 0.4 and xc = 0.6, cr = 0.6. 
Excellent agreement is observed between the final and test mode shapes. Moreover, by comparing 
the final data curve in Figs. 5 and 6 with the mode shape in Fig. 11(c) of Ref. [10], a more accurate 
third mode is approximated. This indicates that more accurate information on mode shapes is 
required to obtain a satisfactory solution from the identification process in the case where both 
crack position and crack depth are unknown. 

Questions arise concerning the conditions under which the identification procedure can pro- 
vided an unique solution. As discussed in Shen and Pierre [10, 11] and concluded in the studies 
of Gladwell et. ai [8], if all the mode information is used in the identification procedure, then 
the system’s properties can be identified uniquely. However, for practical reasons, in structural 
dynamic testing only a small subset of the eigenvalues and eigenvectors can be represented in the 
measurement data. Furthermore, even if substantially more modal information would be avail- 
able, the minimization search may be prohibitive for such a large-dimensional feasible domain that 
would result. These comments are intended to point out certain limitations inherent in the identifi- 
cation procedures. These considerations is addressed with the presentation in the following, which 
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describes sufficient conditions for the unique identification from the dynamic measurements of a 
multi DOF vibrating spring-mass system. 

Spring-mass model 

The following examples of damage identification problems were constructed by introducing the 
damage through the drop in the stiffness or, more conveniently, the value of each damage parameter 
to change the system’s dynamics behaviour. These dynamic changes, taken as the test simulation 
of response data, are used to deduce the value of each damage parameter via the identification 
process. 

The numerical optimization technique set forth in this study for vibrating cracked beam iden- 
tification problems is accomplished using the VMCON optimization package program (this imple- 
ments a sequential quadratic programming method). The damage properties = 1,2,3) of 

the fixed-free spring-mass system are identified by direct solution of the optimalization problems 
described in the previous section. 

The first example corresponds to the identification of a system’s damage, <fmi = 0.0, dm 2 = 
0.5 , dmz = 0.25, using first and second mode information. The first five variables, £i,£ 2 , dm\, 
dm 2 , and dm3 of each vector x are listed in Table 3. The top row denotes the assumed damage 
parameters and corresponding first and second eigenfrequencies and the symbol * denotes the 
expected optimal solution through the identification process. The first and second column entries, 
indicate the fundamental and the second frequencies corresponding to the initial damage 
parameters, dmi , dm 2 , dm3, which are given in the next three columns. The last five columns give 
the final values corresponding to previous entry values. These final values are obtained at the stage 
of the program is terminated when the further optimal search obtain improvements F less than a 
tolerance (10E — 5 was adopted in the present study). 

In Table 3, each case has the final values of £ close to but almost all of them have the 
unacceptable final results for xc and cr. Only the case with initial dm 2 — 0-48 and dm3 = 
0.24 has less than a 1% estimation error. These results show performance of the present damage 
identification process is generally unacceptable if only first and second modes are used. 

The first six variables, (1,(2* £3* dm\, dm 2 , and dm3 of each vector x are listed in Table 4 the 
top row denotes the assumed damage parameters and corresponding first and second eigen frequen- 
cies and the symbol * denotes the expected optimal solution through the identification process. In 
this example, all the modes are used to deduce the damage conditions. Satisfactory predictions are 
obtained in each case, in contrast to the results examined in Table 3. Even though starting point 
is located at boundary of the feasible set (dmi = 0.0, dm 2 = 0.0,d7?i3 = 0.0), the agreement is still 
precise. These results confirm the expectation that a unique and accurate solution predictions are 
assured if all the modal information is included as data in the damage identification process. 

Conclusions 

A general method for damage identification of a simple beam and a spring-mass system is 
presented. The method may be useful as a component of an on-line nonintrusive damage detection 
technique for vibrating structures. A formulation is expressed as a direct minimization problem 
statement with a criteria of the mean square difference of natural frequencies and mode shapes 
between test measurements and corresponding model values. The damage identification problem 
is reduced to finding the damage parameters that will satisfy appropriate constraints and minimize 
the mean square difference. 

The uniqueness and reliability of the identification process is confirmed by solving several 
damage identification examples with specified damage positions. Without knowing the damaged 
location, a restricted region in initial data space had been found for which there will be a realistic 
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if substantially more modal information would be available. However the minimization searc > 
be prohibitive for such a large-dimensional feasible domain that would resu . 
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| Test Data; £7 =0.84703. £7 = 70.1348. cr‘=0.5 | 

initial Data 

Final Data 

1 6 

& | ct 

£1 !'fc t cr | 

1.0 

81.0 i 0.0 

0.84684 | 70.1348 1 0.50033 | 

0.98841 

80.0769 j 0.1 

0.84697 1 70.1346 , 0.50019 i 

0.97217 

PT8.8135 i 0.2 

0.84704 | 70.1347 I 0.49998 1 

0.94815 

77.0062 i 0.3 

0.84701 | 70.1348 1 0.50007 | 

0.91032 

74.3024 | 0.4 

0.84694 | 70.1347 [ 0.50024 j 

0.73638 

i 63.7848 j 0.6 

0.84705 i 70.1348 [ 0.49962 j 

0.54574 

i 55.0511 I 0.7 

0.84703 1 70.1348 1 0.50034 j 

0.27233 

45.9316 f 0-8 

0.84700 j 70.1347 | 0.50009 1 


Table 1: Numerical results based on mean square problem statement of Eqs. (21-25) with the crack 
(damage) specified ( xc = 0.5). 



Test Data: 

=0.84703. Q 

= 70.1348, cr'=0.5 

, ic’=0.5 


Initial Data 

Final Data i 

T 


6 

cr 


ft 

ft 



9 

0.91806 

78.5161 

0.4 

0.4 

0.69639 

70.1359 

0.99789 

0 .36289 1 

9 

0.91371 

76.6365 

0.4 

0.43 

0.70007 

70.1362 

0.99440 

0.39620 j 

9 

0.91158 

75.1335 

0.4 

0.46 

0.84610 

70.1347 

0.91029 

0.53775 i 

9 

0.91056 

74.7464 

0.4 

0.47 

0.84711 

70.1347 

0.67125 

0.49033 j 

9 

0.91063 

74.5157 

0.4 

0.48 

0.84704 

70.1348 

0.50554 

0.49972 | 

9 

0.73472 

63.8062 

0.6 

0.51 

0.84704 

70.1348 

0.60027 

0.50526 | 

9 1 

0.73711 

64.2643 

0.6 

0.52 

0.84704 

70.1348 

0.60083 

0.50531 i 

9 i 

0.73617 

64.7619 

ITT" 

0.53 

0.84704 

70.1348 

0.60141 

0.50534 | 

9 i 

0.73929 

65.6727 

faJ - 

0.54 

0.84705 

70.1348 

0.60255 

0.49459 j 

9 j 

0.73909 

66.6112 

0.6 

0.55 

0.84702 

70.1348 

0.99721 

0.24709 i 

9 ! 

0.75452 

! 74.0109 

0.6 

0.6 

0.70040 

70.1363 I 

0.99079 

| 0.59307 | 

45 1 

0.97475 

80.2193 

032 

0.4 

0.90130 

70.1347 

0.94855 

| 0.94404 | 

45 1 

1 0.91806 

78.5161 

0.4 

0.4 

0.84420 | 

70.1345 1 

0.53053 

I 0.51586 1 

45 | 

0.91531 

77.2676 

[ 0.4 i 

0.42 

0.84686 | 

70.1347 j 

0.50838 

0.50198 I 

45 I 

0.96219 

78.5819 

0.25 

0.45 

0.84643 1 

70.1348 1 

l 0.51729 1 

I 0.49389 | 

45 1 

0.75452 

74.0109 

0.6 

0.6 i 

0.84645 | 

70.1348 j 

0.51723 

0.50609 t 

45 1 

0.64083 | 

77.7173 


0.7 | 

0.89079 ! 

70.1347 | 

0.5SS95 

0.31317 1 


Table 2: Numerical results based on mean square problem statement of Eqs. (21-25). The position 
of the damage xc is a variable. 


| Test Data; (7= 

0.15296. fj» 

1.2956 

£7=2.2494, dm;. 

=0.0. dm;=0.5. dm^=0.25 

i initial Data 

Final Data | 

fi 

1 6 

am^ 

dm- j 

dm? 

(i 

ii 

dm\ 

dm? 

dm? J 

0.19806 

1 1.5549 

0.0 

0.0 

0.0 

0.15299 

1.2956 

0.21392 

0.40077 

0.10985 1 

0.18986 

| 1.4975 

0.0 

0.1 

0.05 

0.15294 

1.2955 

0.17532 

0.41741 

0.132651 

0.18123 

1 1.4429 

0.0 

0.2 

0.1 

0.15296 

1.295G 1 

0.13478 1 

0.43493 

0.15777 ! 

0.17218 

1.3911 

0.0 

0.3 

0.15 

0.15297 

1.2956 

0.09222 

0.45445 

0.18530 j 

0.16275 

j 1.3420 

0.0 

0.4 

0.2 

0.15293 

1.2955 

0.04751 

0.47635 

0.2157C 1 

0.15846 

1 1.3088 

! 0.0 

0.44 

i 0.24 | 

0.15294 

1 1.2956 | 

0.02043 | 

i 0.48975] 

0.23502 | 

0.15494 

i 1.3047 

0.0 

0.48 

1 0.24 i 

0.15295 j 

j 1.2956 1 

i 0.00973 I 

1 0.49511 I 

0.24281 i 


Table 3: Numerical results for spring-mass model using first and second mode information. 
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Tett D ata: ^ = 0 . 15296 . s0, °‘ 

Lrutui Dtu ! 


final Data 


i ami l am-i \ am-^ | 


<2 


am, 


am 7 


jmi 

0.2498 i 


0.19806 I 1.3549 I 


0.18986 I 1.4975 


3.2469 I 0.0 0.0 i 0-0 


‘fl. 15294 T 1.2966 l 3.2494 I 0.0002 J.5U01 


3.0209 | 0.0 ; 0. 


0.05 i 0.15294 I 132955 > 2.2494 | 0.0002 1 0.5001 


0.2498 r| 
0.24421 


0.18123 I 1.44 29 


2 .sn«fi I n.p : 0.2 I 0.1 i 0.16242 I 1.1353 I 3-2W I 0.0067 ' 0.5026 _l 


0.1721S I 1.3911 


3.6095 I 0.0 i 0.3 


0.16 I 0,1526 5 .131950 i 2.2494 I 0.0039 I 0.5015 i 0.2194 
Q. 15280 i 1.2955 * 2.2494 t 0-0019 I Q.50Q8 l [ * ' 


0.16275 I 1.3420 I 2.4232 I 0.0 1 0-4 1 0-2 


0.15848 i 1.3088 I 2.3329 


0.0 


n.44 i Q.24 i 0.15295 I 1.2956 I 2.2494 1 0.0000 1 Q.5002 


0.2499 il 
1 


0.15494 I 1.3046 i 2.2332 


0.0 


0.48 | Q.24 i 0.15292 I 1.2956 | 2.2494 i 0.0004 : 0.5002 ! 0.2496 


Table 4: Numerical results for spnng*mass model using ail three mode inlormation. 



Figure 1 . CJeomelry of a simply supported beam containing a pair of symmetric edge criuks si 
mid span, re = {. 



2. Geometry of a 3 f>OF spring mui beam modrl 
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